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We study the transverse oscillatory modes of nanomechanical silicon nitride strings under high 
tensile stress as a function of geometry and mode index m < 9. Reproducing all observed resonance 
frequencies with classical elastic theory we extract the relevant elastic constants. Based on the 
oscillatory local strain we successfully predict the observed mode-dependent damping with a single 
frequency independent fit parameter. Our model clarifies the role of tensile stress on damping and 
hints at the underlying microscopic mechanisms. 



The resonant motion of nanoelectromechanical sys- 
tems receives a lot of recent attention. Their large fre- 
quencies, low damping i.e. high mechanical quality fac- 
tors, and small masses make them equally important as 
sensors and for fundamental studies . In either 
case, low damping of the resonant motion is very desir- 
able. Despite significant experimental progress [lol. [Tlj . 
a satisfactory understanding of the microscopic causes of 
damping is not yet achieved. Here we present a system- 
atic study of the damping of doubly-clamped resonators 
fabricated out of prestressed silicon nitride leading to 
high mechanical quality factors • Reproducing the ob- 
served mode frequencies applying continuum mechanics, 
we are able to quantitatively model their quality factors 
by assuming that damping is caused by the local strain 
induced by the resonator's displacement. Considering 
various microscopic mechanisms, we conclude that the 
observed damping is most likely dominated by dissipa- 
tion via localized defects uniformly distributed through- 
out the resonator. 

We study the oscillatory response of nanomechanical 
beams fabricated from high stress silicon nitride (SiN). 
A released doubly-clamped beam of such a material is 
therefore under high tensile stress, which leads to high 
mechanical stability [l3] and high mechanical quality fac- 
tors [l3| • Such resonators are therefore widely used in 
recent experiments [1, Q . Our sample material consists 
of a silicon substrate covered with 400 nm thick silicon 
dioxide serving as sacrificial layer and ah — 100 nm 
thick SiN device layer. Using standard electron beam 
lithography and a sequence of reactive ion etch and wet- 
etch steps, we fabricate a series of resonators having 
lengths of 35/npm, n = {1,...,7} and a cross section 
of 100 • 200 nm as displayed in Fig. la and b. Since the 
respective resonance frequency is dominated by the large 
tensile stress [13, , this configuration leads to reso- 
nances of the fundamental modes that are approximately 
equally spaced in frequency. Suitably biased gold elec- 
trodes processed beneath the released SiN strings actu- 
ate the resonators via dielectric gradient forces to per- 
form out-of-plane oscillations, as explained in greater de- 
tail elsewhere [12| . The length and location of the gold 




FIG. 1. Setup and Geometry a SEM-picture of our sample; 
the lengths of the investigated nanomechanical silicon nitride 
strings are 35/n jim, n = {1, . . . , 7}; their widths and heights 
are 200 nm and 100 nm, respectively, b Zoom- in of a: the res- 
onator (highlighted in green) is dielectrically actuated by the 
nearby gold electrodes (yellow); its displacement is recorded 
with an interferometric setup, c Mode profile and absolute 
value of the resulting strain distribution (color-coded) of the 
longest beam's 4th harmonic as calculated by elastic theory. 



electrodes is properly chosen to be able to also excite sev- 
eral higher order modes of the beams. The experiment 
is carried out at room temperature in a vacuum below 
10"'^ mbar to avoid gas friction. 

The displacement is measured using an interferometric 
setup that records the oscillatory component of the re- 
flected light intensity with a photodetector and network 



analyser [1^ [lj| . The measured mechanical response 



around each resonance can be fitted using a Lorentzian 
lineshape as exemplarily seen in the inset of Fig. 2. The 
thereby obtained values for the resonance frequency / 
and quality factor Q for all studied resonators and ob- 
served modes are displayed in Fig. 2 (filled circles). In 
order to reproduce the measured frequency spectrum, 
we apply standard beam theory (see e.g. [l5[). With- 
out damping, the differential equation describing the spa- 
tial dependence of the displacement for a specific mode 
TO of beam n Mn,m[a;] at frequency /„,i„ writes (with 
p — 2800kg/m"^ being the material density [l3|; Ei, ao 
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FIG. 2. Resonance Frequency and Mechanical Qual- 
ity Factor The harmonics of the nanomechanical resonator 
show a Lorentzian response (exemplary in the inset). Fitting 
yields the respective frequency and mechanical quality factor. 
The main figure displays these values for several harmonics 
(same color) of different beams as indicted by the color. To re- 
produce the resonance frequencies, we fit a continuum model 
to the measured frequencies. We thereby retrieve the elas- 
tic constants of our (processed) material, namely the built-in 
stress (Jo — 830 MPa and Young's modulus E = 160 GPa. 
From the displacement-induced, mode-dependent strain dis- 
tribution, we calculate (except for an overall scaling) the me- 
chanical quality factors. Calculated frequencies and quality 
factors are shown as hollow squares, the responses of the dif- 
ferent harmonics of the same string are connected. 

are the (unknown) real Young's modulus and built in 
stress, respectively) 

Y^Eih'^-^Un,m[x]-ao-^Un,„i[x]-p{2Trfn,in)'^Un,m[x] = 

(1) 

Solutions of this equation have to satisfy the bound- 
ary conditions of a doubly-clamped beam (displacement 
and its slope vanish at the supports {un,m[H / {'^n)] — 
(9/9a;)un,m[±V(2n)] = 0, l/n: beam length). These con- 
ditions lead to a transcendental equation that is numeri- 
cally solved to obtain the frequencies /n,m of the different 
modes. 

The results are fitted to excellently reproduce the mea- 
sured frequencies, as seen in Fig. 2 (hollow squares). 
One thereby obtains as fit parameters the elastic con- 
stants of the micro-processed material Ei = 160 GPa, 
do = 830 MPa, in good agreement with previously pub- 
lished measurements [l3| . 

For each harmonic, we now are able to calculate the 
strain distribution within the resonator induced by the 
displacement u[x] and exemplarily shown in Fig. Ic. The 
measured dissipation is closely connected to this induced 
strain e[x,z,t]_= e[a:, z] exp[i27r/i]. As in the case of a 
Zener model [17| we now assume that the displacement- 
induced strain and the accompanying oscillating stress 



a[x,z,t] = a[x, z]e'}q)[i2'!rft] are not perfectly in phase; 
this can be expressed by a Young's modulus E = Ei + 
iE2 having an imaginary part. The relation reads again 
a[x, z] — {El +iE2)e[x, z]. During one cycle of oscillation 
T = 1//, a small volume 6V of length s and cross section 
A thereby dissipates the energy AUsv = As7ri?2e^. The 
total loss is obtained by integrating over the volume of 
the resonator. 

AUn,m = I dVAUsv = 7^E2 f dFen,m[x,z]2 (2) 
Jv Jv 

The strain variation and its accompanying energy 
loss can be separated into contributions arising from 
overall elongation of the beam and its local bending. 
It turns out that here the former is negligible, de- 
spite the fact that the elastic energy is dominated by 
the elongation of the string, as discussed below. To 
very high accuracy we obtain for the dissipated energy 
AUn,m « -K /l2E2wh^ jidx{d'^ / {dxYun,raY ■ A more rig- 
orous derivation can be found in the Supplementary In- 
formation. 

The total energy depends on the spatial mode (through 
en,m, see exemplary Fig. Ic) and therefore strongly dif- 
fers for the various resonances. To obtain the qual- 
ity factor, one has to calculate the stored energy 
e. g. by integrating the kinetic energy ?7n,m = 

(ixyl/9(27r/„^ni)^Wn,m[a^]^- The overall mechanical qual- 
ity factor is Q = 27rC/n,m/A?7n,m- A more detailed deriva- 
tion can be found in the Supplementary Information. 

Assuming that the unknown value of the imaginary 
part E2 of the elastic modulus is independent of res- 
onator length and harmonic mode, we are left with one fit 
parameter E2 to reproduce all measured quality factors 
and find excellent agreement (Fig. 2, hollow squares). We 
therefore successfully model the damping of our nanores- 
onators by postulating a frequency independent mecha- 
nism caused by local strain variation. Allowing E2 to 
depend on frequency, the accordance gets even better, as 
discussed in detail in the Supplement. 

We now discuss the possible implications of our find- 
ings, considering at first the cause of the high quality fac- 
tors in overall pre-stressed resonators and then the com- 
patibility of our model with different microscopic damp- 
ing mechanisms. In a relaxed beam, the elastic energy 
is stored in the flexural deformation and becomes for a 
small test volume Usv = l/^AsEie^. In the framework 
of a Zener model, as employed here, this result is pro- 
portional to the energy loss (see eq.|3]) and thus yields 
a frequency- independent quality factor Q = E1/E2 for 
the unstressed beam. In accordance with this finding, 
Ref. [13] reports a much weaker dependence of quality 
factor on resonance frequency, in strong contrast with 
the behavior of their stressed beams. 

Similar as in the damping model, the total stored elas- 
tic energy in a beam can be very accurately separated 
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into a part connected to the bending and a part coming 
from the overall elongation. The latter is proportional 
to the pre-stress ctq and vanishes for relaxed beams, refer 
to the Supplement for details. Assuming a constant E = 
El +iE2, Fig. 3 displays the calculation of the elastic en- 
ergy and the quality factor for the fundamental mode of 
our longest {I = 35 pm) beam as a function of overall built 
in stress cto . The total elastic energy is increasingly dom- 
inated by the displacement-induced elongation l7eiong = 
l/2ao'wh /; dx{d/{dx)u[x])'^ . In contrast the bending en- 
ergy C/bend = l/2AEiwh^ Jidx{d^/{dx)^u[x])^, which in 
our model is proportional to the energy loss, is found to 
increase much slower with ctq . Thus one expects Q to in- 
crease with (To , a finding already discussed by Schmid and 
Hierold for micromechanical beams [3]. However, their 
model assumes the simplified mode profile of a stretched 
string and can not explain the larger quality factors of 
higher harmonics when compared to a fundamental res- 
onance of same frequency. Including beam stiffness, our 
model can fully explain the dependence of frequency and 
damping on length and mode index, as reflected in Fig. 2. 
It also explains the initially surprising finding fl^ that 
amorphous silicon nitride resonators exhibit high quality 
factors when stretched whereas having Q-factors in the 
relaxed state that refiect the typical magnitude of inter- 
nal friction found to be rather universal in glassy materi- 
als [i^l • More generally we conclude that the increase in 
mechanical quality factors with increasing tensile stress 
is not bound to any specific material. 

Since the resonance frequency is typically easier to 
access in an experiment, we plot the quality factor vs. 
corresponding resonance frequency in Fig. 3b; with both 
numbers being a function of stress. The resulting rela- 
tion of quality factor on resonance frequency is (except 
for very low stress) almost linear; experimental results by 
another group can be seen to agree well with this find- 
ing f2l'| . In addition, we show in the Supplement that 
although the energy loss per oscillation increases with 
applied stress, the linewidth of the mechanical resonance 
decreases. 

We will now consider the intrinsic physical mechanisms 
that could possibly contribute to the observed damp- 
ing. As explained in greater detail in the Supplement, 
we can safely neglect clamping losses [iJJi thermoelastic 
damping 0, [1] and Akhiezer damping [JQ since they all 
predict damping constants significantly smaller than the 
ones observed. 

Instead, we would like to discuss the influence of lo- 
calized defect states. Similar to the Akhiezer effect, it 
is assumed that the energy spectra of defects are modu- 
lated by strain , which thereby drives the occupancy 
out of thermal equilibrium. In Ref. 27| , this effect is cal- 
culated for a broad spectrum of two-level systems. The 
energy difference of two uncoupled levels as well as their 
separating tunnel barrier height are assumed to be uni- 
formly distributed, leading to a broad yet not flat dis- 
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FIG. 3. Elastic Energy and IVIechanical Quality Factor 
of the Beam in Dependence of Stress a The elastic en- 
ergies of the fundamental mode of tlie beam with I = 35iim 
are displayed vs. applied overall stress separated into the con- 
tributions resulting from the overall elongation and the local 
bending. The dashed line marks the strain of the experimen- 
tally studied resonator ctq ~ 830 MPa, there the elongation 
term dominates noticeably, b Quality factor and frequency 
are calculated for varying stress (tq. In order to compare the 
calculation with other published results quality factor and 
stress are displayed vs. resulting frequency. 



tribution of relaxation rates. In the high temperature 
limit the thus derived energy loss per oscillation becomes 
frequency-independent as assumed in our model. In ad- 
dition, published quality factors of relaxed silicon nitride 
nano resonators [19| cooled down to liquid helium tem- 
perature display quality factors that are well within the 
typical range of amorphous bulk materials [l^j . More- 
over, on a different sample chip we measured a set of 
resonators showing quality factors that are uniformly de- 
creased by a factor of approximately 1.4 compared to the 
data presented in Fig. 2. Their response can still be quan- 
titatively modeled using now an increased imaginary part 
of Young's modulus E^. We attribute this reduction to 
a non-optimized RIE-etch step, leading to an increased 
density of defect states. The corresponding data is pre- 
sented in the Supplementary Information. These three 
arguments clearly favor the concept of damping via de- 
fect states as the dominant mechanism. 

We wish to point out that such a model calculation 
based on two-level systems cannot be rigourously applied 
at elevated temperatures, as the concept of two-level 
systems should be replaced by local excitable systems. 
However, it seems plausible that such a system still ex- 
hibits a broad range of relaxation rates, crucial to explain 
frequency-independent damping. In contrast, a mecha- 
nism with discrete relaxation rates will exhibit damp- 
ing maxima whenever the oscillation frequency matches 
the relaxation rate [1, @, [13] • We further notice that in 
our experiments beams with larger widths exhibit slightly 
higher quality factors. This indicates an increased defect 
density near the surface being either intrinsic or caused 
by the micro-fabrication j28:] (RIE etch). For a fixed 
cross-section however, the applicability of our model is 
not affected. 

In conclusion, we systematically studied the transverse 
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mode frequencies and quality factors of prestressed SiN 
nanoscale beams. Implementing continuum theory, we 
reproduce the measured frequencies varying with beam 
length and mode index over an order of magnitude. As- 
suming that damping is caused by local strain variations 
induced by the oscillation, independent of frequency, en- 
ables us to calculate the observed quality factors with a 
single interaction strength as free parameter. We thus 
identify the unusually high quality factors of pre-stressed 
beams as being primarily caused by the increased elastic 
energy rather than a decrease in damping rate. Several 
possible damping mechanisms are discussed; because of 
the observed nearly frequency-independence of the damp- 
ing parameter E2, we attribute the mechanism to inter- 
action of the strain with local defects of not yet identified 
origin. One therefore expects that defect-free resonators 
exhibit even larger quality factors, as recently demon- 
strated for ultra-clean carbon nanotubes [llj . 
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SUPPLEMENTARY INFORMATION 



DAMPING MODEL 

In a Zener model, an oscillating strain e{t) = ^[e[uj] exp[iujt]] and its accompanying stress a[t] — ^[a[uj] exp[iujt]] are 
out-of phase, described by a frequency-dependent, complex elastic modulus cr(a;) = i?[w]e[a;] = + ii?2[i^])e[w]. 

This leads to an energy loss per oscillation in a test volume SV = 6 A ■ Ss of cross-section SA and length Ss. 

AU5v = J dtEAe[t]--^^{se[t]) ^TrSA5sE2e^ (3) 

force , ^ . 

velocity 

We now employ this model for our case, namely a pre-stressed, rectangular beam of length I, width w and height 
h, corresponding here to the x,y,z-direction, respectively. The origin of the coordinate system is centered in the 
beam. The resonator performs oscillations in the z-direction and, as we consider a continuum elastic model, there 
will be no dependence on the y-direction. For a beam of high aspect ratio I 3> h and small oscillation amplitude, 
the displacement of the m-th mode can be approximately written Um[x, y, z] = itm[a::]. During oscillation, a small test 
volume within the beam undergoes oscillating strain em[x, z,t\. 

This strain arises because of the bending of the beam as well as its elongation as it is displaced. The stress caused by 
the overall elongation is quadratic in displacement, therefore it occurs at twice the oscillating frequency. 

1 / d V 



elongation bonding 

-Aq-^u„Ax]) -(1 + %xpM) + z- 



o ( — ) o (1 + 3^[exp[2icji]]) + z— u„Ja;]3fi[exp[ia;i]] (4) 



Inserting this into eq.|3|and integrating over the cross-section w ■ h, the accompanying energy losses can be seen to 
separate into elongation and displacement caused terms. 

wh f d V wh^ ( d"^ \^ 

AC/„.,i = 7rsi;2[2a;]— ( — j ^ tx sE2\^\—^ \-^u„^\x\\ (5) 

Integrating over the length yields the total energy loss of a particular mode AC/ = /_;/2 dx/!^Uwh- In the case that 
E2 is only weakly frequency-dependent, it turns out that for our geometries the elongation term is several orders of 
magnitude (10^ — 10^) smaller than the term arising from the bending. The energy loss therefore may be simplified 
and writes 

At/ « AC/bonding = 7r_E2^^ j |^ ^-j [x] j (6) 

ELASTIC ENERGY OF A PRE-STRESSED BEAM 

A volume 5V subject to a longitudinal pre-stress (Tq stores the energy Ugv when strained; Ei is assumed to be 
frequency independent in the experimental range (5-100 MHz) 

Usv ^ sA(^aoe+]^Eie^^ (7) 

To apply this formula to the case of an oscillating pre-stressed beam, we insert eq.|4|t=o (maximum displacement) and 
integrate over the cross-section to obtain 

1 il f d \^ 1 f d'^ 1 f Q 
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Analog to eq.[5]we can omit the first term in the brackets; integrating over the length yields 




elongation bending 



We can therefore divide the total energy into parts arising from the elongation and the bending of the beam. Depending 
on the magnitude of the pre-stress, either of the two energies can dominate as seen in Fig. 3a of the main text. We 
have checked that the kinetic energy f/kin = l/2/5(wm)^ da;(wm[x])^; (cj„i/(27r), p: resonance frequency, material 
density, respectively) equals the total elastic energy, as expected. 

FREQUENCY-DEPENDENT LOSS MODULUS 

There is no obvious reason that the imaginary part of Young's modulus E2 should be completely frequency- 
independent. We therefore assume that E2 obeys a (weak) power-law and chose the ansatz: 

E2[f]^E2{f/hr (10) 

Fitting our data with the thus extended theory, we achieve a very precise agreement of measured and calculated 
quality factors, as seen in Fig. SI. The resulting exponent is a = 0.075; E2 varies therefore by 20% when / changes 
by one order of magnitude. 




FIG. 4. Resonance frequeneies and quality factors of the resonators a Measured quality factor and resonance frequency of 
several harmonics of beams with different lengths (color-coded) are displayed as filled circles (same data as in Fig. 2 of the main 
text). The resonance frequencies are reproduced by a continuum model; we calculate the quality factors using a model based 
on the strain caused by the displacement. In contrast to Fig. 2 of the main text and Fig. S2 we here allow E2 to be (weakly) 
frequency-dependent. 



LINEWIDTH OF THE MECHANICAL RESONANCE 

The elastic energy of a harmonic oscillator is given hy U = l/2rncffi^o^o with mcff' , , xq being effective mass, 
resonance frequency and displacement, respectively. If we assume the effective mass to be energy- independent, it 
applies Wo oc VU- Recalling the definition of the quality factor Q = 2ttU /AJ7 oc U, one obtains for the for the Full 
Width at Half Max (FWHM) of the resonance 

VU At/ 

As in the main text, the energy depends on the applied overall tensile stress. FigurelH] shows a numerical calculation 
of the resulting linewidth vs. applied stress; one can see that increase in energy loss per oscillation is dominated by 
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the increase in energy, resulting in a decreased linewidth. The exact effective mass is included in this calculation; as 
it changes by less than 20%, the above assumption is justified. 




Prestress [Pa] 



FIG. 5. Linewidth of the mechanical resonance The calculated linewidths (FWHM) for the fundamental mode of the beam 
with / = SSym are displayed vs. applied overall stress. 



MICROSCOPIC DAMPING MECHANISMS 

We start with clamping losses as discussed, e.g., in ref. P, i- e. the radiation of acoustic vifaves into the bulk 
caused by inertial forces exerted by the oscillating beam. With a sound velocity in silicon of vsi ~ 8km/s, the 
wavelength of the acoustic waves radiated at a frequency of 10 MHz from the clamps into the bulk will be greater 
than 500 pm, and thus substantially larger than the length of our resonators. Considering each clamping point as 
a source of an identical wave propagating into the substrate, one would expect that mostly constructive/destructive 
interference would occur for in-/out-of-phase shear forces exerted by the clamping points, respectively. With clamping 
losses being important, one would therefore expect that spatially asymmetric modes with no moving center of mass 
exhibit significant higher quality factors than symmetric ones [3| . Another way to illuminate this difference is that 
symmetric modes give rise to a net force on the substrate, whereas asymmetric modes yield a torque. Since the 
measurement (Fig. 2) does not display such an alternating behavior of the quality factors with mode index m (best 
seen for the longest beam), clamping losses are likely to be of minor importance. 

The next damping mechanism we consider are phonon-assisted losses within the beam. At elevated temperatures, 
at least two effects arise, the first being thermoelastic damping: because of the oscillatory bending, the beam is com- 
pressed and stretched at opposite sides. Since such volume changes are accompanied by work, the local temperature 
in the beam will deviate from the mean. For large aspect ratios as in our case, the most prominent gradient is in the 
z direction. The resulting thermal flow leads to mechanical dissipation. We extend existing model calculations (i| to 
include the tensile stress of our beams. Using relevant macroscopic material parameters such as thermal conductivity, 
expansion coefficient and heat capacity we derive Q- values that are typically three to four orders of magnitudes larger 
than found in the experiment. Therefore, heat flow can be safely neglected as the dominant damping mechanism. 
In addition, the calculated thermal relaxation rate corresponds to approximately 2 GHz, so the experiment is in the 
so-called adiabatic regime. Consequently, one would expect the energy loss to be proportional to the oscillation 
frequency, in contrast to the assumption of a frequency independent E2 and our experimental findings. 

Another local phonon-based damping effect is the Akhiezer-effect Q ; it is a consequence of the fact that phonon 
frequencies are modulated by strain, parameterized by the Griineisen tensor. If different phonon modes (characterized 
by their wave vector and phonon branch) are affected differently, the occupancy of each mode corresponds to a different 
temperature. This imbalance relaxes towards a local equilibrium temperature, giving rise to mechanical damping. 
In a model calculation applying this concept to the oscillatory motion of nanobeams [6| , the authors find in the case 
of large aspect ratios length/height that the thermal heat flow responsible for thermoelastic damping dominates the 
energy loss by the Akhiezer effect. We thus can safely assume this mechanism to be also negligible in our experiment. 
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REDUCED QUALITY FACTOR 

We fabricated a set of resonators, shown in Fig. Sla, that showed lower quality factors than the ones presented in 
the main text (Fig. 2); we attribute this reduction to a non-optimized RIE-etch step. As in the main article, it is 
possible to reproduce the quality factors using a single fit parameter, namely the imaginary part of Young's modulus 
E2- The ratio of the two sets of quality factors is displayed in Fig. SI b and can be seen to be around 1.4 with 
no obvious dependence on resonance frequency, mode index or length. A non-optimized etch step causes additional 
surface roughness and the addition of impurities, thereby increasing the density of defect states. As there is no obvious 
reason why another damping mechanism should be thereby influenced, we interpret this as another strong indication 
that the dominant microscopic damping mechanism is caused by localized defect states. 
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FIG. 6. Comparison of the resonance frequencies and quality factors of the sets of resonators a Measured quality factor and 
resonance frequency of several liarmonics of beams with different lengths (color-coded) are displayed as filled circles. The 
resonance frequencies are reproduced by a continuum model; a model based on the strain caused by the displacement allows 
us to calculate the quality factors, shown as hollow squares. The uniform reduction of the Q-factors is attributed to an non- 
optimized RIE-etch. b The ratio of the quality factors of the two sets resonators (Fig. 2 main article and Fig. S2a) are displayed 
vs. frequency, being approximately constant. 
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